Geometric frustration and magnetization plateaus in quantum spin and 

Bose-Hubbard models on tubes 
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We study XXZ Heisenberg models on frustrated triangular lattices wrapped around a cylinder. 
In addition to having interesting magnetic phases, these models are equivalent to Bose-Hubbard 
models that describe the physical problem of adsorption of noble gases on the surface of carbon 
nanotubes. We find analytical results for the possible magnetization plateau values as a function of 
the wrapping vectors of the cylinder, which in general introduce extra geometric frustration besides 
the one due to the underlying triangular lattice. We show that for particular wrapping vectors 
(A'^, 0), which correspond to the zig-zag nanotubes, there is a macroscopically degenerate ground 
state in the classical Ising limit. The Hilbert space for the degenerate states can be enumerated by 
a mapping first into a path in a square lattice wrapped around a cylinder (a Bratteli diagram), and 
then to free fermions interacting with a single Zjv degree of freedom. From this model we obtain 
the spectrum in the anisotropic Heisenberg limit, showing that it is gapless. The continuum limit 
is a c = 1 conformal field theory with compactification radius R = N set by the physical tube 
radius. This result cannot be checked against a Lieb-Schultz-Mattis argument, for the argument is 
inconclusive when applied to this problem. We show that the compactification radius quantization 
is exact in the projective J±/Jz <^ 1 limit, and that higher order corrections reduce the value 
of R. The particular case of a (A'^ = 2, 0) tube, which corresponds to a 2-leg ladder with cross 
links, is studied separately and shown to be gapped because the fermion mapped problem contains 
superconducting pairing terms. 

PACS: 67.70+n, 71.10.Pm, 75.10.Jm 



I. INTRODUCTION 

There has been growing interest in one-dimensional 
quantum spin chains and ladders that display magneti- 
zation plateaus when subjected to an external magnetic 
field 0-0]. This requires materials with exchange cou- 
plings within a range for which the necessary external 
fields are attainable in the laboratory. Examples are 
the spin 5=1 nickel compound [Ni(Medpt)2(/i-ox)(^- 
N3)]ClO4-0.5H2O, which displays a plateau at (M) = 1/2 
1^, and the spin 5=1/2 compound NH4CUCI3, which 
displays plateaus with magnetization (Af) = 1/4,3/4 ||] 
(the magnetization is measured as a fraction of the max- 
imum magnetization per spin S). One of the most in- 
teresting features in these systems is the possibility of 
gapless plateaus even in integer spin systems |^]. 

An important ingredient to obtain the plateaus is a 
p-merization, or the presence of periodic structures that 
allow for ground states with non-zero (M) other than the 
trivial cases (M) = ±1. A particular case of a p-merized 
structure is a "spin tube" , like the one proposed theoret- 
ically by Cabra, Honecker, and Pujol and by Citro, 
Orignac, Andrei, and co-workers [^j^. These tubes basi- 
cally consist of p-leg spin ladders with periodic boundary 
conditions coupling the 1st and pth chains. Although in- 
teresting structures from a theoretical perspective, these 
tubes are not yet realizable experimentally. 

In this paper we explore another type of spin tube, 
which is formed by wrapping a triangular lattice on a 



cylinder. The triangular lattice with anti-ferromagnetic 
couplings is frustrated, and, depending on the wrapping 
vector, extra geometric frustrations are introduced. One 
motivation to study such types of spin tube lattices is 
that they are realized physically in monolayer adsorp- 
tion of noble gases on the surface of carbon nanotubes. 
This paper contains a detailed account of the results in 
Rcf. |l^, as well as new analytical results that explain 
the numerical findings of the previous work. 




FIG. 1. An example of wrapping of the graphite sheet to 
make a (3, —2) tube. a± are the primitive lattice vectors of 
the honeycomb lattice. The solid rectangle is the supercell of 
the tube. Also shown is the tripartite labeling A, B, C. 

The lattice sites for the spin tubes we consider can 
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be described by a pair of integers {N,M), exactly the 
notation for carbon nanotubes, for which the fohow- 
ing nomenclature applies: {N,N) are called armchairs, 
{N, 0) are zig-zags, and the general case are chiral. The 
triangular lattice points are the centers of the hexago- 
nal lattice of the nanotubes (the adsorption centers), as 
shown in Fig. |l|. The integers {N, M) define the wrap- 
ping vector of the honeycomb lattice that identifies the 
origin with the point Na+ + Ma- = 0, where a± are the 
primitive lattice vectors of the honeycomb lattice. 

The tripartite nature of the triangular lattice of ad- 
sorption sites is destroyed whenever the wrapping vector 
(7V,M) is such that iV - M is not divisible by 3. An 
example of a (7,0) tube is presented in Ref. [0, and 
reproduced below in Fig. 0. 
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FIG. 2. Adsorption sites on a (7, 0) zig-zag nanotube 

The problem of adsorption of noble gases onto planar 
graphite can be understood within a lattice gas model, 
with the hexagonal substrate providing adsorption sites 
that form a triangular lattice . This problem is equiv- 
alent to a Bosc-Hubbard model, where a large nearest 
neighbor repulsion arises from the fact that even the 
smaller of noble gas atoms, helium, spread over an area 
larger than that of the carbon hexagons (the character- 
istic length scale for the zero-point motion of the He 
atom comes from the Lennard- Jones potential length 
scale (7=2. 56A, while the separation between C atoms 
in the hexagons is 1.42A) [Q. The potential energy for 
atoms to be adsorbed in the same site are even larger, 
and thus we take the hard-core boson limit. The case of 
adsorption on the nanotube geometry is similar once the 
periodicity condition due to the wrapping vector (iV, M) 
is considered. 

The lattice gas is defined by the Bose-Hubbard Hamil- 



where rii is the boson density at site i, V is the nearest 
neighbor repulsion and t is the hopping amplitude. In 
the equivalent Heisenberg spin representation, 

n=-J±Y.(^^s'^+^f^')-J^T.^?^^'^T.^^ (2) 

{ij) {ij) i 

where Sf = n^- 1/2, Jj_ = 2t, = -V, and H = ^I-W 
is an effective external magnetic field. The equivalence 
between the hard core Bose-Hubbard model and the XXZ 
Heisenberg spin S* = 1/2 model allows us to study both 
physical problems of the magnetization properties of the 
spin tube and the adsorption of noble gases on carbon 
nanotubes at the same time. 

In the adsorption problem one finds filling fraction 
plateaus as a function of the chemical potential; this 
is the counterpart of the magnetization plateau in the 
spin models as a function of external magnetic field. The 
filling fraction plateaus most often correspond to solid 
phases, with an energy gap in the excitation spectrum. 
However, as we show in this paper, there are examples of 
gapless plateaus which occur in the case of {N, 0) tubes 
(zig-zags) with N not divisible by 3. These plateaus cor- 
respond to a compressible correlated fluid, which is de- 
scribed by a c = 1 conformal field theory with compacti- 
fication radius R ^ N. 

Many of the interesting results we find in the study of 
the spin tubes come from the geometric frustration that 
is introduced when the lattice is wrapped up into a cylin- 
der. In some cases, one can think of the boundary con- 
ditions as forcing a single domain wall running along the 
tube. Recently, Henley and Zhang ||l5| have studied the 
nearest neighbor Hubbard model for spinless fermions as 
a toy model for understandying stripe phases when the 
system is doped away from half-filling. The situation is 
very similar to our case. 

The paper is organized as follows. In section || we 
study the classical limit when the hopping t in the Bose- 
Hubbard model goes to zero; in spin language, this corre- 
sponds to the Ising limit. We explain a macroscopic de- 
generacy present in the frustrated {N, 0) zig-zag tubes. 
In section [H we turn on the quantum hopping term t 
(or the XY spin couplings J± = 2t). We sumarize our 
findings in section |v| where we also discuss possible ex- 
perimental signatures of the plateaus. 



II. CLASSICAL LIMIT 

We start our analysis in the simpler classical limit, 
which corresponds to the Ising limit J± = 0. In the 
Bose-Hubbard model it is the no-hopping limit t = 0. 
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We will carry out the discussion interchangeably between 
the Bose-Hubbard and spin languages. 

Without loss of generality, we take the integers N, M > 
0; all other cases can always be brought to this form by 
a rotation of the basis vectors a± by multiples of 60°. 
The triangular lattice sites can be partitioned into three 
interpenetrating sublattices that we label by A, B, and 
C. If we define q = {N — M)(mod 3), then the wrapping 
of the lattice does not break this tripartite property if 
q = 0. 

Let us begin by studying the commensurate case (g — 
0) before we continue to the more interesting incommen- 
surate cases {q = 1,2). As a function of the chemical po- 
tential, the possible filling fractions are v = 0,1/3,2/3, 
and 1, corresponding to the cases of no filling, one sub- 
lattice filled, two sublattices filled, or all three sublattices 
filled. Notice that there is a threefold ground state de- 
generacy for I' = 1/3 from filling any one of the three 
sublattices; similarly, there is a threefold degeneracy for 
ly = 2/3 from filling all but one of the three sublattices. 
The only values of chemical potential for which there is 
a macroscopic degeneracy in the q — case is right at 
the transition points between plateaus. In the spin lan- 
guage, the filling fractions correspond to magnetizations 
(Af) = —1,-1/3,1/3, and 1. A schematic plot of the 
magnetization as a function of the external field H is 
shown in Fig. || 
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FIG. 3. Schematic representation of the magnetization (or 
filling) plateaus for a g = commensurate tube (T = 0). 

Let us now analyze the incommensurate cases q = 1,2. 
One still obtains the trivial plateaus at = 0, 1, but the 
question is what happens to the partially filled plateaus. 
Here we describe an analytical argument that leads to 
the possible filling fractions for general {N,M). In this 
derivation, we make some assumptions based on physical 
intuition, which is supported by exact numerical transfer 
matrix calculations for a range of values for N, M . First, 
one can separate the lattice sites into three classes lo- 
cally, but not globally; there will always be a topological 
line defect running along the tube for q = 1,2. We call 
this line the "zipper" ; examples are shown in Fig. ^ for a 



(4, 0) tube. One can open up the tube along the zipper, 
and label the hexagons centered at the triangular lattice 
points A, _B,and C . However, now the number of A, B, 
and C points are no longer all equal. Filling only one 
of these three classes of points leads to filling fractions 
V that depend on the zipper geometry, and on which of 
A, B, or C were chosen. For example, a case leading to 
filling fractions v — 7/20,3/10 is shown in Fig. |^a, and 
another leading to filling fractions = 1/2, 1/4 is shown 
in Fig. |b. 




FIG. 4. Two examples of zippers for a (4, 0) tube. When 
the tube is opened up along the zipper (thick lines) , the hexag- 
onal lattice can be separated locally into three sublattices 
A, B, and C. By filling the sublattices labelled by A, B, and 
C one obtains for case (a) either u — 7/20 if ^ or B is chosen, 
or !^ = 3/10 if C is chosen. For case (b), one obtains i/ = 1/2 
if A is chosen, or = 1/4 if B or C is chosen. 

There is an infinite number of possible zippers for a 
given {N,M). In order to choose the one that leads to 
the lowest energy state, we argue that the zipper must 
be as straight as possible. The rationale is that bending 
the zipper should cost energy. So zippers as in Fig. ^ 
should be preferred. For a given wrapping vector {N, M), 
N, M > 0, we then investigate two zippers that we label 
by two vectors zi = (1,-2) and Z2 = (2,-1). These 
two vectors differ by a rotation of 60°, and are shown 
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in Fig H A third = (1,1) is not considered because 
for N, M > this zipper winds rapidly around the tube, 
and hence there is a large energy cost associated with its 
longer length. 




FIG. 5. Two straight zippers in the (2, 1) chiral tube, la- 
beled by the vectors (solid arrows) zi = (1, —2) in (a) and 
Z2 = (2, —1) in (b). The wrapping vectors is the dotted ar- 
row. 

To proceed, we must calculate the boundary energy 
cost associated with the topological mismatch on the two 
sides of the zipper. Consider the portion of a tube shown 
in Fig ^ for a (2, —1) zipper, where in each hexagon we 
write its coordinates in the a± basis. 




FIG. 6. Labelling of the hexagons by their coordinates in 
the a± basis for the (2, —1) zipper. The zipper is shown as a 
dark line. The region below the line lies within the boundaries 
of the unwrapped tube. 

Within the boundaries set by the zipper, we can label 
the sites A, B, and C according to: 



(mod 3)^ A 

1 (mod 3) ^ S 

2 (mod 3) ^ C 



(3) 



In Fig. we show the possible boundary configurations 
for q — 0,1,2. Let ni, be the number of near neighbors 
that are now both occupied per unit length, \z2\, along 
the zipper. (We also refer to the situation of occupied 
near neighbor sites as a bond between the sites). For 
q = 0, Ub = since the tripartite lattice is not frustrated, 
but 

• for g = 1, nf, = if -B or C is the filled (minor- 
ity) sublattice; nt = 2 if ^ is the filled (majority) 
sublattice. 

• for g = 2, rifc = if C is the filled (minority) sub- 
lattice; Tifc = 1 if A or S is the filled (majority) 
sublattice. 

Each set of occupied neighbors costs energy V. The en- 
ergy E(v) per area A for a given filling configuration is 



Eiiy) _ ribV 



A A 

where A is the supercell area given by 

A = {Na+ + Ma-) ■ Z2 

= {Na+ + Ma-) ■ (2a+ - a_) 
- i2M + TV) Acx 



(4) 



(5) 



and y^hcx = |a±P'\/3/2 is the area of one hexagon, which 
we will set to unity. 




(a) 





(b) 



(0 



FIG. 7. The boundary configurations for the zipper 
zi — (2,-1) and the cases (a) q — 0, [h) q — 1, and (c) 
q — 2. The hexagons below the solid line lie within the bound- 
aries set by the zipper and are labeled according to eqn. 
Above the solid line, (A^, M) = (0, 0) is labeled by A, while 
{N-1,M + 1), (iV + 1,M) = (1, 0) are labeled by B. 
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Given that the number of A, B, and C lattice points 
are not all equal for q ^ 0, the possible filling fractions 
depend on the choice of filling either majority or minority 
lattice points. It is simple to check that 



A 

[A/31 
A 



for majority filling 
for minority filling 



(6) 
(7) 



These fillings can be converted to magnetization values 



by m± = (M) 



1 and = yu - 3^. Notice 



that the spin flip symmetry m — > —m corresponds to the 
particle-hole symmetry ly ^ 1 — i^. 

From this point on, the analysis is straightforward; one 
has to consider the energy density /A for all the dif- 
ferent possibilities q = 1,2, zippers zi^2, as well as ma- 
jority or minority fillings h'±. The transition between 
and occurs at a critical field He when — E{v^). 

After considering the eight cases, we obtain the min- 
imum energy configurations for a generic (iV, M) tube 
with q= \: 



I 



1 



TO_ = - 1 

2M + N 3 1 2A^ + M 



2M 

N + M 



V 



(8) 



The complementary case of g = 2 is obtained by inter- 
changing N ^ M . We have verified the magnetizations 
and critical fields of eqn. (||) for small values of N, M 
previously [Tc| ]. 

Notice that we recover the graphite (planar) result 
m± = ±1/3 as the tube radius approaches infinity 
(iV, M oo), as we should. 



A. Macroscopic Degeneracy 

Most plateaus for general {N, M) tubes have only a 
discrete degeneracy corresponding to filling different sub- 
lattices. There is a macroscopic degeneracy only at the 
transition points. However, a careful consideration of the 
plateaus in the exceptional case of the frustrated zig-zag 
[N, 0) tubes reveals that either the m+ or m_ plateau is 
macroscopically degenerate for its whole range of stabil- 
ity as a function of chemical potential (or magnetic field). 
Which plateau is degenerate depends on the the value of 

Let us reconsider the filling of the (5,0) tube in 
Fig. ^. Imagine building a typical state layer-by-layer 
from top to bottom, with a total of L layers. Each new 
layer must add exactly two filled sites and one nearest- 
neighbor pair of occupied sites (1 bond). There are two 
ways to achieve this constraint. A non-trivial one is such 
that no two adjacent sites may be occupied within a layer, 
and this is the case shown in Fig. || (left). Notice that 




moving a particle from a layer to another adds another 
intra-layer bond. There is also a trivial one where two ad- 
jacent sites are occupied within a layer, and to conserve 
n;,, two adjacent sites must be occupied in the next, and 
so on down the tube. However, the class of such trivial 
states is only 5-fold degenerate, as opposed to the macro- 
scopically degenerate class containing the state shown in 
Fig.g We now turn to the enumeration of the non-trivial 
class of states. 



(4,2) (2,5) (5,3) (3,1) (1,4) 

\/\/\/\/\ 

(1,3) (3,5) (5,2) (2,4) (4,1) 

/\ / \/\ /\/ 
(4,2) (2,5) (5,3) (3,1) (1,4) 
\ /\/\/ \/\ 
(1,3) (3,5) (5,2) (2,4) (4,1) 

/\ /\/\/\/ 
(4,2) (2,5) (5,3) (3,1) (1,4) 

\/ \/\/\/ \ 
(1,3) (3,5) (5,2) (2,4) (4,1) 

/ \/\/\ /\/ 
(4,2) (2,5) (5,3) (3,1) (1,4) 



FIG. 8. LEFT: Typical configuration in v+ (or m_) of the 
(5, 0) tube. Alternating numbering within layers allows a 
symmetric description from bottom-to-top or top-to-bottom. 
CENTER: Allowed states as paths on a wrapped square lat- 
tice. The vertex labels may be dropped. RIGHT: Steps to 
the right can be represented by particles (solid circles), and 
to the left by holes (open circles). 

An allowed state can be represented as a string of occu- 
pied sites, <T — {ai\, i — I, . . . , L, which in our example 
is (T = {• • • (5, 3)(5, 2)(5, 3) • • •}. At each layer, there are 
exactly two possibilities for the following one. For exam- 
ple, (1, 4) can be followed by (I, 4) or by (2, 4). However, 
the total number of possibilities at any given level is five. 
Fig. 1^ (center) summarizes this structure succinctly as a 
square lattice wrapped on the cylinder. A typical state, 
then, is a lattice path along the tube. A path on this 
type of lattice is known as a Bratteli diagram. There is a 
recent Hubbard model considered by Henley and Zhang 
Ip^ of spinless fermions on a square lattice in which the 
bookkeeping of states is very similar. 

Generalizing to (A^, 0), we find possible states in 
each layer and two in the succeeding one, and the struc- 
ture of states is again that of a wrapped Bratteli dia- 
gram with A^ sites along the circumference. The dimen- 
sion of the Hilbert space is the number of lattice paths, 
A^2^, so that in an infinitely long tube, the entropy per 
site (in the thermodynamic L oo limit) is exactly 
S = (In2)/A^. Notice that constrained paths introduce 
correlations along the length of the tube, despite the ab- 
sence of inter-layer hopping. 

The preceeding discussion is valid for all A^ with q = 2, 
where the macroscopically degenerate plateau occurs at 
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v+. For 5 = 1, the only modification is that the macro- 
scopically degenerate plateau has filling v- . 



III. QUANTUM LIMIT 

Having understood in detail the classical limit, we now 
turn on a small hopping, t V ^ that introduces quan- 
tum fluctuations. Away from the transitions, adding a 
particle necessarily increases n;, so that all plateaus be- 
gin with a classical gap of order V . Consequently, we 
can project to the Hilbert space of the classical ground 
states. Those plateaus which have only a discrete sym- 
metry must retain their gaps, but the macroscopically 
degenerate plateaus are more complicated. 

The matrix elements of the projected Hamiltonian (|l|) 
connect only those states that differ by a single hop of a 
boson: 



{a'\H\a) 



if T.^5a 

otherwise. 



= L-l 



(9) 



where ct, a' are two paths in the Brattcli lattice con- 



structed in Section II A 



It turns out that this Hamiltonian is exactly solvable, 
being closely related to a class of solid-on-solid models 
that were introduced by Pasquier In the following, 
we derive the continuum limit of 7i directly. 

In the previous section, the paths cr were labeled, for 
clarity, by a string of occupied sites on the nanotube. 
A simpler representation is to work with the wrapped 
square lattice directly, where the path is uniquely spec- 
ified by an initial point and its direction in each layer. 
We will label the topmost layer by i = 1 with i increas- 
ing by 1 with each downward move. There are L layers of 
hexagons and we impose periodic boundary conditions, 
L + 1= 1. In order for the layers to match, L has to be 
even. 

To specify the initial point on cr, we chose an "anchor" 
a on one of the N sites in the i — 1 layer (a = 1, . . . , N). 
Now, we represent a step to the right in layer i by a 
fermion, and a step to the left by a hole. We define 
fermion creation and anihilation operators and Ci for 
each layer i. A state jcr) in the Hilbert space, S", is rep- 
resented by 



|0), 



(10) 



where a is the anchor site in the first layer and ii ■ ■ ■ ip 
are the layers where the path steps to the right. For 
instance, the portion of the path in Fig. ^ is jcr) — 
1 3) (8) cjcgcj ••• |0). The fermionic representation is con- 
venient since there is exactly one step in each layer, but 
hard-core bosons can also be used, as in one dimension 
they are equivalent. The number of particles (steps to 
the right) and holes (steps to the left) must add up to L 



(the total number of steps). Each path also has a topo- 
logical character for the number of times that it winds 
around the tube, which must be a multiple of N in order 
for the path to close with the periodic boundary condi- 
tion L + 1 = 1. These two requirements are summarized 
by 



Np 



L 

bN 



(11) 
(12) 



where Np^h is the number of particles or holes, and b is the 
winding number (positive or negative). If each particle 
(hole) is assigned a positive (negative) unit charge, then 
bN is the total charge. Notice that 6 = at half-filling 
{Np ^Nu^ L/2). 

Whenever it is allowed within a layer, a single hop 
changes the step sequence right-left to left-right and vice 
versa, which corresponds to cl_^_^Ci or cjci+i. In a layer 
without a kink no hops are possible, and the hopping 
terms vanish by fermionic statistics. Since we are work- 
ing in periodic boundary conditions, the boundary terms, 
cjci and c|^ci, must be treated more carefully. A hop at 
this point is necessarily accompanied by a translation of 
the anchor point by \a) i-^ |a ± 1). Let us represent this 
operation by 



R±\a) = |a± 1) 



(13) 



with R[_ — i?+ . Cylindrical wrapping requires a Z^r sym- 
metry because \aztN) = \a), i.e. — R±. Putting 
the bulk and boundary hopping terms together, the 
Hamiltonian of eqn. (0) becomes 



n = -~2t 



■i+l'- 



' c\cL 



(14) 



We can think of H as describing free fermions on a pe- 
riodic one dimensional chain with a Z^r impurity on one 
of the bonds. 

7i can be diagonalized exactly in momentum space. 
Going around the tube lengthwise contributes a phase 
gi/cL .^yjjiig going around the perimeter contributes e*"^, 
with (/) = 27Ta/N (a= 1,...,N Therefore toroidal 

boundary conditions require 



Or, 



2ttv 



(15) 



(16) 



where v — 2t \s the velocity and n is an integer. The 
Hamiltonian contains the usual free particle dispersion. 



n 



- cos kc\ck , 

k 



(17) 
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but the allowed k are given by eqn. (|l^). If the spectrum 
(near half-filling) is linearized around the Fermi momen- 
tum, Ikp], then at small k states with nonzero a cost an 
additional energy of {2'kv / L){a/NY . 

The a /N offset in k can be thought of as a minimally 
coupled vector potential such that the magnetic field 
is a ^-flux tube through the torus containing a/N flux 
quanta. In other words, one of the bonds along the chain 
(the anchor) had Z^r symmetry, whose effect is equivalent 
to a flux tube. Fig. || illustrates this equivalence. 



A a ^ 




FIG. 9. The left ring shows the Zjv impurity on the anchor 
bond (wavy line). The right ring shows the equivalent alter- 
native, where the impurity is replaced by a flux tube through 
the torus. 

The offset in k due to the flux induces a total current 
in the fcrmion system. 

At this point, one can see two topological effects on 
the torus. The first is the Zat flux tube, or total cur- 
rent. As we have seen, its contribution to the energy 
near \kp\ was {2TTv/L){a/Nf. The second is the path 
winding along the length of the tube, or the total charge 
defined by eqn. (|l2|). Its contribution to the energy near 
|fci?| is similar, {2-kv / L)[bN /2)'^ . The total energy due to 
these topological sectors is 



2'KV 



a 



(18) 



This expression is familiar from the Luttingcr liquid 
model of one-dimensional spinless Fermions [|l9| . 

We can now obtain the continuum limit of our model. 
It is well known that free fermions in one dimension are 
equivalent to free bosons. The corresponding Lagrangian 
is 



an 



(19) 



where ip is the bosonic field. £ is a conformally invariant 
theory with central charge c = 1. We conjecture that the 
topological effects that we described above come from 
compactifying on a circle of radius R, 



2ttR 



(20) 



By compactifying the boson, topological modes (or zero 
modes) appear. In field theory, they are conventionally 



obtained from electric and magnetic monopoles. The en- 
ergy of the zero modes is 



2ttv 



a 



(21) 



where a and b are integers labeling the fundamental cy- 
cles on the torus. Comparing ( ^ to AEa.b (|T1), 
we find that R = N. The ordinary phonon (oscillator) 
modes exist on top of each topological sector and simply 
contribute the usual phonon energy, so that the complete 
dispersion is 



E 



2tt. 



(22) 



We have verified this energy spectrum by exact diagonal- 
ization for small systems in an earlier work | [lO| . 

The overall picture of a compactified boson with cen- 
tral charge c = 1 is consistent with the solid-on-solid 
models of Pasquier ll^ . One can also ask what happens 
at higher order in t/V, which is the subject of the next 
subsection. 



A. Finite t/V 

The next corrections to H are of order t'^/V and in- 
volve virtual transitions to adatom configurations that 
are not in the degenerate subspace S. The generic form 
is 



— Vs 

V ^ 



(tj)(ki) 



Vs 



(23) 



where Vs is a projection operator into S. Another way 
of writing the second order perturbation is the familiar 
form. 



{a'\-H\\){m\a) 

E\ — Ea 



(24) 



where jcr), \ a') G S while |A) ^ 5 is the intermediate vir- 
tual state with energy E\. E„ and E^i are, of course, 
equal and constant, and all energy differences are due to 
the nearest neighbor repulsion VniUj. 

There are three types of virtual processes: (i) single 
particle hopping from a to a ^ a' , (ii) two particle cor- 
related hopping from a to a ^ a' and (iii) single particle 
diagonal hopping from a back into a. For concreteness, 
consider process (iii) in the (5, 0) state of fig. ^. When- 
ever there is a kink in tr, such as in the third layer from 
the top, the contribution to eqn. (24) from all virtual 
hops is — (35/6)4t^/T^. For example, the adatom on site 
3 can hop into any one of its six neighbors with the en- 
ergy denominators 1/2-^1/2 + 1/2 + 1/2+1/2 + 1/3 
(in units of /V). Similarly, the adatom on site 5 con- 
tributes 1 + 1 + 1, for a total of 35/6 (it is forbidden to 



7 



hop one site over to the right because the resulting state 
is in S). On the other hand, if there is no kink, the con- 
tribution is — 8 • At"^ /V. The criterion for a kink in layer 

i is 2[l/4— {hi — l/2)(ni+i — 1/2)] = 1, where hi — c\ci 

I 



is the fermion number; otherwise this quantity vanishes. 
Similarly, the absence of a kink is synonymous with the 
nonvanishing of 2[l/4 + {h^ - l/2)(ni+i - 1/2)]. Thus, 
the total diagonal contribution to TL can be written 



n->i-L 



V 



EE 



35 
~6 



13 _ 



13 



(25) 



For general N, the correction scales like N. The es- 
sential term in the last line of eqn. (^5|) is the first one. 
This four-fermion interaction rcnormalizes the radius R 
by corrections of order t/V. 

Let us return to processes (i) and (ii). An example of 
(i) is the adatom in the third layer from the top, site 5, 
hopping to the second layer, site 1, and then back into 
the third layer, site 1. The intermediate state is not in 
S. This process serves only to renormalize t because its 
amplitude is the same for all kinks. An example of (ii) is 
the particle in the fourth layer, site 5, hopping to site 4, 
followed by the particle in the third layer, site 5, hopping 
to site 1 in the same layer. This correlated hopping oc- 
curs in a configuration containing the sequence particle- 
hole-hole or hole-particle-particle, which corresponds to 
a next-nearest neighbor interaction cl_^_2cl_^^Ci^iCi + h.c. 
This type of term has a similar effect to the hihi^i term 
in eqn. (p5|), it reduces the value of R by corrections of 
order t/V. The four-fermion terms do not cancel unless 
there is a delicate cancellation. 




FIG. 10. An example of the virtual process (iii) in the con- 
figuration of Fig. 1^ 

One can also consider the extreme limit in which 



I 

t ^ V. In this case, the XXZ Hamiltonian in eqn. (||) is 
simply the XY model on a cylinder. Note that our model 
has the particle-hole symmetry t ^ —t dX least to order 
/V . Let us denote the spin angle relative to the cylin- 
drical surface by 'p(x, 0), where x is the coordinate along 
the tube and 9 is the coordinate around the perimeter. 
Uniqueness of the wavefunction requires that has the 
periodicity ^p{x, 9 -\- 2n) = ip{x, 6) -\- 2TTm, where m is 
an integer. The low energy excitations are purely along 
the length of the tube; excitations around the perimeter 
will cost an energy on the order of 1/A^, which is large 
compared to 1/L. Thus we can freeze the 9 coordinate, 
and the energy density is proportional to IQj^vJp. Since 
the periodicity is still tp = ip + 27rm, we end up with a 
free boson compactified on radius R^^ — 1. The ques- 
tion is how the adsorption regime t <^V, which is also a 
compactified boson but on radius R = N , is reached. 

B. The Lieb-Schultz-Mattis Argument 

Another indication that our plateaus are nontrivial 
comes from trying to understand the spin tube as a spin 
ladder. Let us take this point of view to see if it yields 
the zero gap. A standard approach is to use a Lieb- 
Schultz-Mattis (LSM) argument, in which the spins are 
deformed slowly along the length ^ . Applying it to our 
tube, we find that a plateau is gapless if S — M is not an 
integer, where S and M are the total spin and magneti- 
zation, respectively, per layer (a layer being the N sites 
around the perimeter). Using S = N/2 and the magne- 
tizations from Eqn. (^ , we find that S — M is an integer 
in the macroscopically degenerate plateaus, so that the 
LSM argument is insufficient in this case. Therefore, a 
conclusive argument must take the geometric frustration 
into account. 

C. Special Case: N = 2 

Before concluding with the effective theory, we should 
point out that the geometry of the (2, 0) tube is special; 
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all sites in adjacent layers are interconnected. As a re- 
sult, all of its plateaus have an extensive entropy, and we 
find that hopping opens a gap in both plateaus. Fig. |ll| 
illustrates this exception. At either filling, the adatom in 
each layer is free to hop to either site — both configura- 
tions are iso-energetic because each site is contiguous to 
all sites in the neighboring layers. Hence both plateaus 
are macroscopically degenerate. In the presence of hop- 
ping, each adatom lives in a double well potential, which 
has a gap of order t. Therefore, neither plateau is gapless 
unlike higher N . 

In fact, this tube can be written as a spin chain that 
has been studied at isotropic coupling [Q, —2t = V. 
Two plateaus were found in this case, and it is tempt- 
ing to speculate whether the two regimes are connected 
adiabatically. 



(c) 





FIG. 11. The fillings for iV = 2. u+ = 1/2 (a) and = 1/4 

(b) . Each adatom can live at either site in its layer because 
each site is connected to every site in the neighboring layers. 

(c) shows the triangular lattice in the plane; the horizontal 
double bond is due to the periodicity around a cylinder, (c) 
is exactly the geometry of the spin ladder studied by other 
authors (albeit with different coupling). 

Another way to understand the gap is within the 
fermion model. Due to the high connectivity of the 
Bratteli diagram for ~ 2, there is another allowed 
hop in addition to the c\cij^i -f cj+i'^* term, which is 
c|4-i-i "I" Cj+iCj. This is a pairing interaction that leads 
to a superconducting gap. 



IV. DISCUSSION AND CONCLUSION 

In this paper we studied spin and Bose-Hubbard mod- 
els in tube geometries. We found that an undclying frus- 
trated triangular lattice, combined with extra geometric 
frustration from the closed topology of the tube, leads 
to interesting plateau structures for filling fractions y (in 
the Bose-Hubbard models) or magnetizations (M) (for 
the spin models). In section || we studied the classical 
limit when the hopping t in the Bose-Hubbard model goes 
to zero; in the spin language, this correspond to the Ising 



limit. We obtained the different filling fractions allowed 
energetically by the wrapping condition, labelled by the 
vector ( A^, M) . We showed that the corresponding solid 
phases are unique states, with the exception of the zig-zag 
tubes (A^, 0) with A^ not divisible by 3. In this particular 
case, which we studied in section HA, we showed that 



there is a macroscopic degeneracy of the classical ground 
state for a range of chemical potentials (or fields), i.e. a 
degenerate plateau. In contrast, when A^ = (mod 3), 
there is only a single value of the chemical potential (or 
magnetic field, H = 0) for which the ground state degen- 
eracy is non-trivial. We enumerated all the degenerate 
states at the plateau by mapping each state in the de- 
generate Hilbert space into a path in a Bratteli diagram 
wrapped on a cylinder. Using this enumeration, we then 
turned on the quantum hopping term t (or the XY spin 
couplings J± = 2t) in section 



III. We found that the 



quantum terms lift up the degeneracy, leaving a gapless 
spectrum described by a c = 1 conformal field theory 
with compactification radius R — N . We presented an 
analytical argument that shows that the Bratteli path 
can be mapped into a state of a fermion model in a one 
dimensional chain, and the hopping terms in the tube lat- 
tice correspond to a hopping in the fermion chain. The 
compactification radius follows from a boundary Zjv de- 
gree of freedom and path winding, which must be in- 
cluded when the system is subject to periodic bound- 
ary conditions. We also discussed the cases when the 
anisotropy t/V (or J±/Jz) is no longer small, and the 
special case (A = 2, 0). 

There are a number of interesting features that emerge 
from this problem of restricting spin and Bose-Hubbard 
Hamiltonians to tube geometries. For example, it is 
noteworthy that in the case of the (A, 0) zig-zag tubes 
with A^ not divisble by 3, the low energy spectrum in 
the quantum case is gapless and described by a confor- 
mal field theory with a quantized compactification radius 
i? = A; in the language of Luttinger liquids, this means 
that the Luttinger parameter is fixed by topology, simi- 
larly to the case of edge states in a fractional quantum 
Hall fiuid pO| , and in contrast to quantum wires (where 
the Luttinger parameter can vary continuously). In the 
spin tube problem, this quantization, in the projective 
limit J_L <C \Jz\ (or t <^ V), follows simply from the 
fact that there is a direct connection between the phys- 
ical lattice and the target space for the bosons. One 
of the interesting features of the spectrum for the gap- 
less plateaus is that the specific heat is linear in tem- 
perature T, with the prefactor related to the velocity 
and a universal constant that depends on the central 
charge c = 1: C = c^^^T = ^^^^ T. This dependence 
would be manifest in thermal measurements, and could 
be contrasted with exponentially activated behavior for 
the gapped solid phases. 

Recently, there has been some interest in a connec- 
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tion that was pointed out by Kitaev between the 
topological stability in some class of spin Hamiltonians 
in toroidal structures and quantum error correction for 
quantum computation. Bonesteel p2| investigated the 
possibility of using the two-fold topological degeneracy 
of spin- 1/2 chiral spin hquid states on the torus to con- 
struct quantum error correcting codes. Since closed car- 
bon nanotube structures (tori) have been observed ex- 
perimentaly [ p3[ , it may be possible that spin tube ge- 
ometries such as the ones discussed in this paper could 
realize physically some related Hamiltonians. We note, 
however, that for the class of Hamiltonians with nearest 
neighbor couplings that we studied here, there was only a 
trivial non-degenerate ground state, although there were 
topologically non-trivial excitations. 

Finally, we would like to discuss a possible experimen- 
tal set up for observing the filling fraction plateaus in the 
case of the monolayer adsorption problem of rare gases 
on the surface of carbon nanotubes. The filling fraction 
plateaus could be probed by measuring the resonance fre- 
quencies for a vibrating single wall carbon nanotube as a 
function of the vapor pressure of rare gas in a chamber. 
One way of measuring the resonant frequencies would 
be to have a single-wall tube hang alongside a multi- 
wall tube, both with metallic grains at their tips. The 
multiwall tube is very stiff, and is basically rigid as com- 
pared to the single wall. By measuring the changes in 
the capacitance between the grains, one could probe the 
frequency of vibrations of the single- wall tube. This fre- 
quency is given by the square root of the ratio between 
the linear mass density and the stifness of the tube. Dif- 
ferent filling fractions of adsorbed noble gases change the 
mass density of the tubes; however, their effect on the 
tube stiffness should be insignificant. The change in mass 
density of the nanotube due to adsorption of helium is 
basically 5p/ p — v M-ii(./2Mc = v/12. Hence, the ratio 
between the resonant frequencies W1.2 for two plateaus at 
1/12 is 



(26) 



Resonant frequencies can be measured very precisely, so 
the filling fraction steps vi — V2 should be experimentally 
measurable as a function of the vapor pressure of helium. 
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